Tuning the atomic and domain structure of epitaxial films of multiferroic BiFeC>3 
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Recent works have shown that the domain walls of room-temperature multiferroic BiFeOs (BFO) 
thin films can display distinct and promising functionalities. It is thus important to understand 
the mechanisms underlying domain formation in these films. High-resolution x-ray diffraction and 
piezo-force microscopy, combined with first-principles simulations, have allowed us to characterize 
both the atomic and domain structure of BFO films grown under compressive strain on (001)- 
SrTiOs, as a function of thickness. We derive a twining model that describes the experimental 
observations and explains why the 71° domain walls are the ones commonly observed in these films. 
This understanding provides us with a new degree of freedom to control the structure and, thus, 
the properties of BiFe03 thin films. 

PACS numbers: 



Magnetoelectric multiferroics exhibit coupled electric 
and magnetic orders, which might lead to a variety of 
novel devices that would benefit from the fact that the 
magnetization (polarization) of these materials can be 
controlled by means of an electric (magnetic) field [l|. 
For practical devices, multiferroics are preferred in thin 
film form. Moreover, the strain induced by the mismatch 
between the film and the substrate lattice parameters 
can sometimes be used to tune the film properties with 
respect to the bulk Q. 

Bismuth Ferrite, BiFeOs (BFO), is one of the few mul- 
tiferroics that orders above 300 K and, thus, one of the 
most promising ones The ferroelectric properties of 
BFO are very robust, and it displays record polarization 
values of about 100 /iC/cm 2 . Since the ground state of 
bulk BFO is rhombohedral (space group R3c), symmetry 
arguments suggest that the thin films grown on cubic sub- 
strates under compressive epitaxial strain should be mon- 
oclinic (space group Cm or Cc, depending on whether the 
06 octahedra rotations are clamped by the substrate or 
not, respectively). Indeed, several authors 0, @, 0] 
have reported a monoclinic unit cell that is similar to 
that of strong piezoelectric PbZri-^Ti^Oa (PZT) with 
x « 0.5 Q. The proposed link between the strong piezo- 
electricity and the symmetry of the unit cell [9j], which 
allows the polarization to rotate, adds to the interest of 
BFO films 
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Beyond their intrinsic properties, BFO films are cur- 
rently receiving renewed attention because of the novel 
functionalities observed to occur at domain walls (DWs). 
Indeed, recent works have shown that some BFO DWs 
are highly conductive [10[ , and that the DW density con- 
trols the magnitude of the (exchange bias) coupling be- 
tween BFO and other (metallic) layers in complex het- 
erostructures [11]. It is thus of prime importance to 
achieve control of the domain structures and understand 
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their formation. In contrast, it is striking to note the 
scarcity, and lack of agreement, of experimental informa- 
tion on the atomic structure of the films and its evolution 
with thickness 0, @, [l2[. Indeed, we believe that a com- 
plete picture of the structure of these films does not exist 
yet. 

We have grown BFO thin films on SrRu03-buffered 
SrTiOs (STO) substrates, and followed the unit cell dis- 
tortion as a function of thickness during the first stages 
of strain relaxation. Our c/a ratios are consistent with 
those in Refj6|- Additionally, we have been able to resolve 
the monoclinic distortion and measure the evolution of 
the full unit cell. The comparison of the experimental 
results with several structural models simulated ab initio 
allowed us to resolve the monoclinic space group (Cc) and 
atomic structure, as well as the polarization direction. 

Several (OOl)-oriented BFO thin films with thickness 
ranging from 12 to 87 nm were grown on atomically fiat, 
Ti02-terminated (001)-STO substrates with low miscut 
angle (0.1°). Conductive layers of SrRuOa with a thick- 
nesses of 5 nm were deposited in between the substrate 
and the BFO layer. The BFO films were grown by pulsed 
laser deposition, assisted by reflective high energy elec- 
tron diffraction (RHEED), using a KrF excimer laser 
(A= 248 nm). The deposition was performed at 670° C 
in an oxygen pressure of 0.3 mbar. After deposition, the 
films were cooled down slowly to room temperature un- 
der an oxygen pressure of 100 mbar. 

The evolution of the crystallographic distortion with 
thickness was investigated by mapping the reciprocal 
space using x-ray diffraction (XRD) from lab (out-of- 
plane) and synchrotron (in-plane) sources. Due to the 
epitaxy, which fixes the [001] direction in reciprocal space 
to be perpendicular to the substrate surface, the twelve 
possible monoclinic domains are reduced to four and the 
reciprocal space maps are significantly simplified (similar 
to the case of a crystal under an electric field [13]). In 
particular, if all domains are present, looking at the ar- 
eas around the substrate (hhl) c , which corresponds to the 
(h01) m reflection in the monoclinic structure [lj^ one can 
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FIG. 1: a) Monoclinic domains in BFO thin films under com- 
pressive strain, in the (h01) m (left) and (Okl) m (right) scatter- 
ing planes, b) Diffraction map around the (h01) m reflections 
when all four domains are present, c) Directions of the polar- 
ization and the monoclinic lattice parameter, a m , for the four 
down polarized domains 



extract the three lattice parameters and the monoclinic 
angle, as sketched in Fig. 1. 

Typical reciprocal space maps (RSMs) around the 
(113) c STO substrate reflections for ultrathin (< 18 nm) 
and thin (> 18 nm) films are shown on Fig. 2a. The 
RSMs of the thinnest BFO films display a broad (113) c 
peak (using the pseudo-cubic notation), at the same Ky 
(in-plane component of the scattering vector) of the sub- 
strate, showing that they are fully coherent with the 
substrate. The FWHM of these films agrees with what 
is expected for their thicknesses. There is thus no in- 
dication of unresolved splitting. The RSMs of thicker 
films display a splitting of the (113) pc BFO peak, as ex- 
pected (see Fig. 1). The monoclinic lattice parameters 
extracted from these patterns are plotted in Fig. 3. Inter- 
estingly, c m shows no changes with increasing thickness. 
This is in agreement with the report by Kim et al. [6[, 
who showed that the lattice parameters of the strained 
films are constant below ~ 100 nm, a puzzling and un- 
explained result. However, for thickness above 18 nm we 
observe a splitting of the in-plane parameter values and a 
/^90°, characteristic of a monoclinic distortion. Figure 3 
reveals a gradual increase of the monoclinic distortion 
&m — bm with thickness. In addition, grazing incidence 
XRD (GID) has shown that the in-plane pseudo-cubic 
angle, j pc is, indeed, different from the out-of-plane an- 
gle (3, and that such a difference decreases with increas- 
ing thickness (7 pc = /3 in the relaxed structure). Inter- 
estingly, the deviation of a m and b m from the value of 
2xd( 110 ) (i.e., the fully coherent case) is symmetric. 

First-principles simulations [15[ allowed us to ratify 
these results and gain further insight into the atomic 
structure of the BFO films. It has been shown that fer- 
roelectric thin films can be successfully studied by sim- 
ulating the corresponding bulk material subject to elas- 
tic boundary conditions that mimic the epitaxial con- 



straints imposed by the substrate. In this work we ex- 
tend such an approach to make a distinction between the 
cases of ultra-thin and thin films, for which we consider 
different elastic constraints. More precisely, in the ultra- 
thin (uth) case we assumed the film is strongly clamped 



by the substrate, and impose a p 



asTO and 



7 pc =90°. In contrast, in the thin (th) case 



we only imposed that the in-plane area be constrained 
to be a| TQ . This allowed us to model the ultra-thin to 
thin transition evidenced by the experimental results of 
Fig. 3. On the other hand, in our simulations we consid- 
ered two structural models, with and without rotations 
of the 06 octahedra, which correspond, respectively, to 
the Cc and Cm space groups. 

Our simulations clearly indicate that the BFO films 
present significant 06 rotations and thus the Cc space 
group. Indeed, when allowing for 06 rotations we com- 
puted 4V«pc h=1 - 0016 and ^c/^c h=0 - 9981 for the split- 
ting of in-plane lattice parameters, in good agreement 
with the experimental values of 1.0015 and 0.9980 de- 
rived from the data in Fig. 3. In contrast, when the 
06 tiltings are clamped in the simulations, we obtained 
4c/ a pc h =i- 0029 and ^/^c h=0 - 9968 - The c pc /a pc ra- 
tios follow the same pattern: the value computed with 
(without) 06 tiltings is about 1.07 (1.17), to be com- 
pared with the experimental result of approximately 1.04, 
which strongly suggests that even in the thinnest films 
the 06 rotations are not fully clamped by the substrate. 
These results provide a justification to first-principles 
studies of monoclinic BFO films in which a structural 
model with 06 rotations is adopted (see, e.g., Ref. [liH ). 
Additionally, for the calculated monoclinic angle we ob- 
tained /3=90.4°, which is perfectly compatible with our 
experimental results, and we computed c^/c^ h = 1.0005, 
in agreement with our experimental observation that the 
c pc lattice constant is weakly dependent on thickness. Fi- 
nally, the computed polarization is very weakly affected 
by the uth-to-th transition: We obtained P= 90 /iC/cm 2 , 
with in-plane and out-of-plane components of 59 jaC / cm 2 
and 68 /iC/cm 2 , respectively. The polarization forms an 
angle of about 11.6° with the body diagonal of the pseu- 
docubic cell, being rotated towards the [001] direction. 

Let us now describe the evolution of the domain struc- 
ture. Analysis of the RSMs in Fig. 2 shows that the 
domain walls that prevail are the 71° ones (see Fig. lc) 
in agreement with previous reports [17[. This is con- 
firmed by grazing incidence diffraction (not shown here) , 
as well as from PFM images (Fig. 2b). Out-of-plane 
PFM measurements show that all the films are polar- 
ized down, also in agreement with previous reports [I?}. 
In Fig. 2b, in-plane PFM (IP-PFM) images of the same 
films are shown. In agreement with the XRD data, we 
observe a clear evolution of the domain pattern. For the 
thinnest films, no contrast is detected on the IP-PFM 
images. IP-PFM images for intermediate films show a 
clear stripe-like pattern. These stripes indicate four po- 
larization variants, which are in good a gree ment with 
rhombohedral-like monoclinic distortions [18|](see figure 
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FIG. 2: a) Reciprocal space maps around the (103) m (= (113) pc ) reflections for different thicknesses of the BFO films. The 
axes are in units of 2K = 47r/A, with A= 1.5405A. The horizontal line through the maps indicates the out-of-plane reciprocal 
lattice spacing, which is unchanged in the range of studied thicknesses, b) In-plane piezo-response images of the same films. 
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FIG. 3: Evolution of monoclinic lattice parameters and shear 
displacement as a function of thickness. 



lc) and with the RSM maps. We observe that the num- 
ber of variants decreases from four to two variants with 
further increasing thickness, allowing for longer stripes 
for the thicker films. Several works have already shown 
two- variant stripe domains for (QQl)-oriented BFO films, 
by using miscut STO substrates [12|, [l9[ or orthorhombic 



substrates [20j. The origin of this reduction of polariza- 
tion orientations in BFO films was reported to be the 
step-flow growth and the substrate anisotropy, respec- 
tively. Since the growth mode as well as the substrate 
miscut in all our films are the same, our results point to 
yet a different mechanism. 

All this evidence fits a simple but powerful model by 
which the domain formation enables and controls the 
monoclinic distortion of the unit cell. Figure 4a shows 
how twinning reduces the in-plane strain introduced by 
the pseudo-cubic angle, j pc (characteristic of the mon- 
oclinic distortion). Two pairs of twins, coherent along 
[100] (vi) or along [010] (V2), can form. It can also be 
seen that, in order to do that, the in-plane lattice pa- 
rameters of the film, df and 6/, deviate equally from 
the fully strained values of a srt = b str = V^sto, i.e., 
Sf = a str — S and bf = b str — as sketched in Fig. 4b. 
The magnitude of the shear vector 5, therefore, deter- 
mines both cif and 6/, which split symmetrically with in- 
creasing thickness (see Fig. 3). As a result, the in-plane 
area of the film is unchanged with respect to the fully 
coherent film, which in turn seems compatible with the 
observation that the c m lattice parameter does not vary 
during strain relaxation (for thicknesses up to 100 nm). 
A more subtle result of this relaxation process is that the 
symmetry of film unit cell is actually lower than mono- 
clinic; indeed, it can be seen from Fig. 4b that the angle 
between df and bf is given by 7/ = cos _1 (£ 2 /(ifbf), and 
thus different from 90°. A very similar twining mecha- 
nism with symmetry lowering has been found in thin films 
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FIG. 4: a) The two types of twins present in the films, each 
including 71° walls: vi (left) is coherent along [100] and V2 
(right) is coherent along [010]. b) Detail of the film distortion. 



of TbMn0 3 grown on (001)-STO substrates [2l[, which 
suggests it may be typical of low symmetry perovskites on 
cubic substrates. As observed in figure 4a), the two pairs 
of variants, 90° rotated from each other, are in agree- 
ment with the PFM maps of the 42 and 59 nm films and 
give rise to both 71° and 109° walls. Even though both 
nucleate with equal probability in the growing film, on 
a low miscut substrate, because of the relatively large 
strain energy store at the boundary between them, for 
thicker films (and therefore larger strain energy at those 



boundaries) one of the two variants will be preferred, as 
observed in the thicker 89 nm film. In the presence of a 
substrate miscut, the steps can indeed determine which 
of the two variant is present [12], but this will happen 
in exact substrates provided that the films have enough 
time to relax. 

In summary, we have observed clear trends in the evo- 
lution with thickness of the structure and microstructure 
of BiFe0 3 films on (001)-SrTiO 3 substrates. We have 
shown that the lattice parameters and the film symme- 
try do not result simply from the mismatch with the sub- 
strate, but also from the occurrence of a particular twin- 
ning that allows for the observed monoclinic distortion. 
This twinning model explains why the 71° domain walls 
are so often observed in atomically fiat films on (001)- 
SrTiOs. Such an effect provides us with a new degree 
of freedom for tuning the structural and physical proper- 
ties of the thin films. Our results suggest that the physics 
behind the effects of epitaxial strain is richer than usu- 
ally thought, and that traditional thermodynamic phase 
diagrams and first-principles models need to be comple- 
mented with knowledge of the domain structure in order 
to reach a full understanding of the materials behavior. 
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